Spin coherence lifetime extension in Tm^+:YAG through dynamical decoupling 
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We report on spin coherence lifetime extension on Tm^"^:YAG obtained through dynamicaUy 
decouphng the thulium spins from their magnetic environment. The coherence lifetime reached with 
a Carr-Purcell-Meiboom-Gill sequence revealed a 450-fold extension ['-^(230 ± 30) ms] with respect 
to previously measured values. Comparison to a simple theoretical model allowed us to estimate the 
correlation time of the fluctuations of the ground level transition frequency to (172±30) at 1.7 K. 
For attaining eflicient decoupling sequences, we developed a strategy inspired by the zero-first-order 
Zeeman effect to minimize the large inhomogeneous broadening of the ground level spin transition. 
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I. INTRODUCTION 

The building block of quantum information is the su- 
perposition state between two states of a quantum sys- 
tem, known as a qubit. The success of a certain qubit 
technology depends strongly on how long the quantum 
character of the state holds. Reaching long coherence life- 
times of the quantum system is thus a main issue. This 
is even a more crucial concern when it comes to quantum 
memories, whose goal is to keep the information as long 
as possible, respecting faithfully the quantum character 
of the stored qubit f^i^ Rare-earth- based systems are ex- 
cellent prospects of quantum memories thanks to their 
long coherence lifetime, T2, at low temperature, both 
of optical and hyperfine transitions. In particular, rare- 
earth ion doped crystals (REIDC) present the additional 
advantage of providing motion-less centers without any 
trade-off loss of T2.^ Therefore, rare-earth long T2S can 
be fully exploited in REIDC since the storage time of the 
device is not limited by atomic motion as it is in atomic 
vapors. Even though storage protocols make use of the 
longer-lived ground-state hyperfine (spin) superposition 
states, extension of T2 from the ~100-/Lts to the ^100-ms 
scale, or even the second, is desirable in REIDC. 

At low temperature, decoherence is ruled by the inter- 
action of the rare-earth spins with the fluctuating mag- 
netic field produced by the nuclear spins of the host ma- 
trix. Different techniques have been developed to decou- 
ple the system from its environment in the seek of extend- 
ing T2. The one known as spin-locking consists in apply- 
ing a strong radio-frequency (rf) magnetic field parallel 
to the spins and tuned to the spin transition. "^'^ However, 
good tuning to the spin transition is not available in in- 
homogeneously broadened systems, such as REIDC. Be- 
sides, when the spin coherences result from the transfor- 
mation of optical ones -as in quantum storage protocols- 
the initial parallelism between spins and rf field cannot 
be ensured for all the spins in the ensemble. Therefore, 
spin-locking extension of T2 is not compatible with quan- 
tum storage. 

A successful technique recently implemented^'^ on 
REIDCs makes use of the zero- first- order Zeeman (ZE- 
FOZ) effect. It consists in carefully looking for a static 
magnetic field precisely sized and oriented such that the 



gradient of the ground level splitting with respect to the 
magnetic field coordinates vanishes. As a result, the spins 
are insensitive (in first order) to the environment- induced 
magnetic field fluctuations. This has led to T2 extensions 
of up to 3 orders of magnitude/^, Nevertheless, ZEFOZ 
critical points are available only in REIDCs that present 
a zero-field hyperfine splittingi^ 

Dynamical decoupling (DD) is a broadly used tech- 
nique for the extension of T2 (for a review, see Ref. Q)- 
Based on the Hahn spin echo;^ it was first applied in 
the context of NMR for high precision spectroscopy 
More recently, the interest for preserving the coherence 
in quantum information qubits has been pointed out»^^— 
The simplest, and still highly efficient, DD scheme is the 
Carr-Purcell-Meiboom-Gill sequenceJ^ii^ It consists in 
the repetition of an elementary sequence composed by a 
T/2-long free evolution interval, a 7r-pulse and another 
T/2-long free evolution interval [see Fig. [SJa)]. The 
phase of the 7r-pulses alternates between and tt in order 
to compensate for pulse imperfections. If the spin level 
splitting does not change significantly during the dura- 
tion T of the elementary sequence, the spin coherence will 
be found at known states at the end of each elementary 
sequence, i.e. sX t = jr, with j the elementary sequence 
index. These states are either the initial state, for j even, 
or a TT-phase-shifted state, for j odd. Hence, the coher- 
ence can be preserved as long as the 7r-pulses are applied. 
The challenge consists in increasing the vr-pulse rate in 
order to compensate for the effect of the environment 
before the latter suffers from significant reconfiguration. 

Experimental realizations of DD have yielded remark- 
able results, especially in moleculeSfii trapped ionsiS 
and nitrogen- vacancy centers in diamond Fraval and 
coworkers have successfully applied the CPMG sequence 
to Pr3+:Y2Si05^ the most thoroughly studied REIDC 
for quantum memory applications. A combination of the 
ZEFOZ and CPMG techniques yielded a coherence life- 
time over 30-s long. 

Thuhum-doped YAG (Tm^+iYAG) has been identified 
as an attractive REIDC thanks mainly to its diode-laser- 
accessible optical line (793 nm) and to its very simple 
level structure under static magnetic field (two ground 
and two excited sub-levels). Its large optical inhomoge- 
neous broadening (20 GHz) can also turn out useful for 
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broadband applications such as light storagej2£ In recent 
years, it has been actively investigated in the prospect 
of quantum storage^"— However, the degeneracy of its 
hyperfine sublevels due to quenching of the total angular 
momentum^^ does not allow the ZEFOZ technique. T2 
extension in Tm'^+ : YAG must thus rely on DD only. Nev- 
ertheless, the simple CPMG sequence cannot be straight- 
forwardly transposed to Tm^+:YAG. The reason, as we 
will discuss in detail below, is its large spin transition 
inhomogeneous broadening (^20 times broader than in 
Pj.3+.Y2Si05), which imposes high power constraints to 
the rf field. Hence, a strategy for minimizing the spin 
inhomogeneous linewidth is compulsory. 

This article reports on T2 extension through CPMG- 
based DD on Tm^+:YAG. In Sec. [Hi we will discuss 
the problem related to inhomogeneous broadening and 
we will describe a ZEFOZ-inspired strategy for minimiz- 
ing the broadening. Three experimental sections follow. 
In Sec. nil Al we will characterize the system (resulting 
inhomogeneous broadening, non-extended T2, etc). In 
Sec. nil Bl we will present results on the viability of the 
CPMG sequence on Tm^+iYAG. Finally, Sec. HITCl win 
be dedicated to T2 extension through DD. The Appendix 
describes a theoretical model that fits our experimental 
data. 



II. MINIMIZATION OF THE 
INHOMOGENEOUS BROADENING 

A CPMG sequence is based on the application of suc- 
cessive rf TT-pulses to a spin ensemble. In order to act 
over the whole ensemble, the hard pulse spectrum must 
cover that of the spins. Thus, rather short pulses are 
necessary. At the same time, the condition of pulse area 
equal to tt must be fulfilled. Therefore, the broader the 
atomic spectrum, the higher the requirements on pulse 
intensity. 

As we said, the CPMG sequence has been successfully 
applied to Pr'^+:Y2Si05. In that crystal, the inhomoge- 
neous hyperfine linewidth, T^nh^ is as narrow as 15-30 
kHz.'^"''^^ A square 7r-pulse of Rabi angular frequency Q, 
has a bandwidth of the order of O. In the case of Tinh 
of the order of some tens of kHz, reasonable rf powers 
can produce rf-pulses able to act on a bandwidth sev- 
eral times broader than Tinh- Pulses need not be shorter 
than a few tens of microseconds. In other words, efiicient 
TT-pulses are readily affordable for Pr'^+:Y2Si05. 

The scenario is quite different for Tm'^+iYAG. As re- 
ported in previous experiments;^ Fin/i ~ 500 kHz, which 
makes rf TT-pulses much more challenging. Indeed, for the 
experimental conditions of Ref. [2^, the 280- kHz Rabi fre- 
quency obtained with a 65-W rf power yields a TT-pulse 
bandwidth far too narrow to expect an efficient CPMG 
sequence. Hence, a strategy for minimizing Tinh is to be 
developed for Tm^+iYAG. 

Let us first address the reason behind this large 
Tinh- Splitting of hyperfine sublevels is achievable in 



Tm^"*" : YAG through the application of an external mag- 
netic field B. For a given B, the energy difference be- 
tween sublevels is determined by the gyromagnetic ten- 
sor 7. It is actually the strong anisotropy of this tensor 
in Tm'^+iYAG (71, ^ 1x,1y', x, y, z to be defined later) 
that is responsible for high Tinh - Let us describe in a few 
words how this happens (we will go through the details 
later): the high anisotropy of the gyromagnetic tensor en- 
tails a strong sensitivity of the sublevel splitting angular 
frequency A to the orientation of B. Either spatial in- 
homogeneities in the field or variations in the lattice ori- 
entation within the sample make that the magnetic field 
orientation seen by the Tm'^+ ions varies slightly from 
one ion site to another. Because of the high anisotropy 
of the gyromagnetic tensor, these slight variations entail 
strong ones in A. This results in a span of values for 
A that add up to give a line broadening. The strategy 
we propose to minimize this broadening is to seek for a 
particular orientation of B where the sensitivity of A to 
the orientation is minimum. 

Although driven by a different motivation, this ap- 
proach is quite similar to the ZEFOZ technique. In that 
case, the procedure consists in searching for a specific B 
whose orientation and magnitude satisfy VsA = O^ili^ 
where Vb is the gradient with respect to the magnetic 
field coordinates. In the problem we are interested in, 
of course, an increase of T2 through the ZEFOZ effect 
would be welcome, but our main concern is still to mini- 
mize Tinh in order to maximize the efficiency of TT-pulses. 
The procedure here also searches for a minimization of 
the sensitivity of A to B. However, full vanishing of V^A 
is possible only in systems that present a zero-field hyper- 
fine splitting.^ This is not the case of Tm'^+iYAG, where 
hyperfine sublevels are degenerate because of quenching 
of the total angular momentumj^ Thus, a full ZEFOZ 
effect is not possible in Tm'^+iYAG. Nevertheless, partial 
ZEFOZ, that is, vanishing derivatives of A with respect 
to two of the magnetic field coordinates, is still possible, 
as we will see in detail below. 

In the presence of a magnetic field, the Hamiltonian 
for Tm^"*" ions takes a simple form for sites of D2 point 
symmetry, which is the case of Tm^"*" substituting for 
Y'^+ in a YAG matrix. It can be expressed as a nuclear 
Zeeman term^ 



n' = -hjnBeff ■ I, (1) 

with an effective magnetic field given by 



Bejy = {ixB^yi + lyByy + 7^S^z)/7„. (2) 

The quantities I and 7„ are the nuclear spin operator 
and gyromagnetic ratio, respectively, and x, y, z are the 
three orthogonal twofold axis of the D2 point group, 
7j,, 7^ are defined as the eigenvalues of the rank 2 tensor 
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7a/3 = 7«1 + 
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where is Bohr's magneton, gj is the Lande factor for 
the manifold of total angular momentum J, Aj is the 
hyperfine interaction parameter and En is the energy 
of the eigenstate \n) of the Hamiltonian comprising the 
free ion and crystal field terms. In the case of D2 point 
group symmetry, 7^^ is diagonal in the basis {x, y, z}i24 
The splitting between the eigenstates of the Hamiltonian 
(H)) is given by 
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FIG. 1. (Color online). Orientations of the six orientationally 
distinct sites of the Y"^"*" ion in the YAG crystal latticei^ Each 
parallelepiped represents the local D2 symmetry for one of the 
sites. The x, y, z axis are the local ajcis for site 4. 



where B = {B^ + By + i?^)^/^, and 9 and (j) are spherical 
coordinates for the polar and azimuthal angles, respec- 
tively, with respect to some convenient Cartesian refer- 
ence frame. 

ions in the YAG lattice occupy six sites that are 
equivalent in terms of their local environment but dis- 
tinct in terms of the orientation of this environment, i.e. 
the axis of the D2 point group are differently oriented 
(see Ref. [35j and references therein). The site orienta- 
tions are sketched in Fig. [T] by six parallelepiped that 
represent the possible orientations of the D2 symmetry. 
When an external field is applied, its effect on each site 
must be treated separately. To study the behavior of 
Ab(0, (/)), we will consider the site for which x, y, z are 
oriented along [Oil], [Oil], [100], respectively, this is the 
site "4" as labeled in Fig. [T] In Fig. [DJa) the strong 
variations of Ab as a function of {6, </>) reflect the high 
anisotropy of 7^, 7y, jz- For the figure, we used the val- 
ues -jy = (403 ± 3) MHz/T, 7^ = 0.0457^, 7^ = 0.0177^ 
compatible with the experimental observations of Ref. 
[21I . The angles 6 and have been taken with respect to 
the axis {[100], [DIG], [001]} of the YAG Bravais lattice. 

Our aim is to minimize Tinh- For that, our strategy is 
to look for a magnetic field orientation {9, (f) for which 
the splitting sensitivity to the magnetic field is minimum. 
We see in Fig. [UJa) that Ab(0, 0) displays both maxima 
and minima that meet this criterion. The question now 
is what exact field orientation to choose. For example, is 
it better to choose a maximum or a minimum? Figure 
[llja) suggests the maximum is a better option since the 
second derivative is lower (in absolute value) than that 
at the minima. However, that is not the only criterion to 
take into account. 

As we said before, the aim of keeping Tinh small is to 
lower the constraint on 7r-pulse bandwidth, and thus on 
TT-pulse duration, so that the pulses produced with the 
available rf power will act efficiently over the whole Tinh ■ 
Of course, rf-power is not the good parameter to look at; 
it is rather its associated Rabi frequency D, which comes 
into play. Usually, it is the rf field component orthogonal 
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FIG. 2. (Color online), (a) Spin sublevel splitting per T 
and (b) Rabi frequency per mT as a function of the ori- 
entation of the static magnetic field B with respect to the 
{[100], [010], [001]} reference hame of the YAG crystal. The 
rf field Bi is always orthogonal to B, and contained in the 
(001) plane. The lines indicate the orientation chosen for the 
experiments: (6 = 54.8°, </> — 45°), which ensures low Tinh 
and high Q. 



to the static field that determines il. For Hamiltonian 
terms such as (H)), both for B and for the rf field Bi, 
what counts is the component of ^i,eff orthogonal to 
Beff (Bi,eff is defined as in Eq. ^ by substituting B^, 
By, Bz by Bi x, Bi^y, i3i,z). In our experimental setup, 
the fields B and Bi are generated by either coils per- 
pendicular to one another. As the coils are fixed, the 
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static field orientation is actually performed by rotating 
the crystal rather than the magnetic field. Therefore, 
the orientations of B and Bi cannot be varied indepen- 
dently. When choosing the orientation of B (or rather 
that of the crystal) care must be taken since, because of 
the strong gyromagnetic tensor anisotropy, D, may vary 
significantly. 

In Fig. IHb) we present the behavior of 51 as a function 
of the angles 9 and (j) that determine the orientation of B 
and, consequently, that of Bi. The experimental setup 
fixes the rf field orientation as Bi/Si = (— sin (p, cos (p, 0) 
with respect to the crystal frame {[100], [010], [001]}. It is 
clear from the figure that a non-negligible f2 is observed 
only at very specific orientations. The region of high fl 
coincides with the "C" -shaped minimum line observed 
for Ab{9, 4>) in Fig. [D^a). This is of course not a coinci- 
dence, it results from the particular anisotropy oijab- To 
understand this statement we need to bear in mind that, 
as soon as a given field has a non- negligible y component, 
its associated effective field is (almost) parallel to y be- 
cause of 7y » 7a;, 7z [see Eq. ^]. Thus, even though 
B ± Bi, it happens that 'B^ff \\ Bi^ejf for most field ori- 
entations. f2 is determined by the component of Bi 
orthogonal to ^eff- The orthogonality only happens if 
either (?) i3y ~ or (ii) Bi y ~ 0. In case (ii), however, 
[Bi ejyl is strongly diminished as a result of ^x,lz '^ly 
Therefore, only case (i) displays high Q.. At the same 
time, as 7y ^ lx,lz, fulfilling (z) results in minimizing 
A [see Eq. That is why the regions of minimum A 

and high coincide in Fig. [2] 

The arguments just given answer the question on 
whether to choose a maximum or a minimum of Ab(0, <j)) 
when aiming at minimizing Tinh- To further specify the 
orientation (0, (j>) of B, it is useful to recall that, at 
(j) = 45°, site 6 presents the same A and as site 4 
described above. This permits to double the experimen- 
tal signal by addressing two sites at a time, and, at the 
same time, to avoid a double-valued fi. Once (/) — 45° 
set, 9 is automatically fixed by selecting the minimum 
of Ab(6',45°). It occurs at 6* = 54.8°. For the result- 
ing orientation [0 = 54.8°, = 45°), indicated in Fig. 
[5] by crossed lines, we expect As = 15.3 MHz/T and 
n/{2TTBi) = 101 kHz/niT. Static and rf fields of the order 
of a few T and mT, respectively, are perfectly compatible 
with the experimental capabilities. From previous exper- 
iments performed at [6 = 49°, (j) = 45°), we estimated the 
angular uncertainty in our crystal to be ~0.3°. The cor- 
responding Tinh obtained theoretically is 28 kHz. The set 
of values we have found for Tinh , A and make CPMG 
sequences experimentally feasible on Tm'^+iYAG. 

As regards the remaining sites, they will not partici- 
pate in the experiment. Sites 1, 3 and 5 display a spin 
transition frequency far from resonance for the chosen 
field orientation. In addition, sites 3 and 5 present a 
forbidden optical dipole moment for the electric field po- 
larization (parallel to [111]) we will use for initializing 
and reading out the system. That is also the case for site 
2. 



III. EXPERIMENTS AND DISCUSSION 
A. Characterization 

The experiments were carried out in a 0.1% at 
Tm'^+iYAG crystal cooled down to 1.7 K in a liquid He- 
lium cryostat. Superconductive coils generate a static 
magnetic field which lifts the spin degeneracy, resulting 
in a four-level system composed of two ground states and 
two excited states. The static magnetic field is oriented 
as described in the previous section. As in Ref [28l. the 
spin transition is resonantly driven by a phase-controlled 
rf magnetic field. A 10-turn coil, 20-mm long and 10- 
mm in diameter conveys the magnetic excitation to the 
crystal sitting at its center. The rf signal generated 
by an arbitrary wave generator (Tektronix AWG5004) 
passes through a pulsed amplifier (TOMCO BT00500- 
AlphaSA). A ~600-kHz bandwidth rf resonant circuit 
further amplifies the rf current fed to the coil. The static 
magnetic field magnitude, close to 1 T, is chosen to pro- 
duce a ground level spitting tuned to the rf circuit res- 
onance (15.1 MHz). The light beam, emerging from an 
extended cavity diode laser and time-shaped by acousto- 
optic modulators, is resonant with the optical transition 
of the system at 793 nm. It is used for initializing and 
probing the system. It propagates along the [110] di- 
rection with polarization parallel to [111]. The optical 
length of the L = 5-mm-long sample has been measured 
to be aL = 0.9. 

First of all, we will experimentally verify the reduction 
of Tinh and the availability of a reasonably high D,. A 
hole burning experiment [see Fig. [2Ia)] revealed a full 
width at half maximum of the spin transition antihole 
of 105 kHz. This is well below the previously reported 
value of c±500 kHz,^* showing that the strategy for mini- 
mizing Tinh has been fruitful. In the absence of inhomo- 
geneous broadening, the antihole width should coincide 
with that of the central hole. The latter is measured to 
be 93 kHz. The ^10-kHz growth from the hole to the 
antihole is consistent with the expected value of Tinh- 
The hole width by far exceeds the homogeneous width of 
the optical transition, of the order of a few kHz only.— 
Since the hole has been deeply burnt, down to the trans- 
parency level, saturation may explain the large observed 
broadening, a typical feature in similar situations. The 
laser linewidth, reduced to less than 1 kHz over 10 ms 
by locking to a Fabry-Perot interferometer, should not 
affect the hole width, but the drift of the reference cav- 
ity during the burning process might also add to hole 
and antihole broadening. An imperfect orientation of B, 
differing slightly from {9 — 54.8°, (/) = 45°), might also 
contribute to Tinh larger than expected. In any case, our 
experimental result fixes an upper bound for Tinh ■ As re- 
gards fi, we have performed a spin nutation experiment 
[see Fig. M.'b)] that yielded n/(27r) = 264 kHz for a rf 
power (^27 W) well below the maximum capability of 
our amplifier. With this value, square-shaped 7r-pulses 
as short as 1.8 ^s are possible. This gives a pulse band- 
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FIG. 3. (a) Hole burning spectrum. Ground level splitting 
antiholes (Ag), excited level splitting holes (Ae) and mixed 
antiholes (Ag — A^) are observed in addition to the central 
hole. The Ag antihole is 105-kHz wide, (b) Spin nutation 
experiment. The rf field is turned on an f = 0. The spins 
nutate at a Rabi angular frequency of 27r x 264 kHz. 



width, '^f2, comfortably larger than Ti„h- These are the 
parameters to be used in the sequences of the CPMG 
experiments to be discussed. 

In order to weigh the impact of the CPMG sequence 
on the extension of T2, we have measured the latter by 
means of an optically detected spin echo experiment. For 
the optical detection, a 7r/2-pulse is applied at the end 
of the spin echo sequence (see Fig. [SJa) for a sketch; 
consider n = 1). It is aimed at converting the coherence 
Pab back into an optically detectable population differ- 
ence. The transmission of a probe beam is measured 
prior and post application of the rf spin echo sequence. 
The relevant experimental information is contained in the 
comparison between the initial and final probe intensity. 
The results are presented in Fig. S) We observe a non- 
linear behavior, the non-linear character being more pro- 
nounced at short times. This tells us that, at the time 
scale of Fig. SI we have not yet reached the limit t ^ Tc 
in which \pab\ decays in an exponential fashion (tc the 
correlation time of the fluctuations in A induced by the 
environmental aluminium spins). We also note that |pq6| 
does not reach unity for times approaching to zero. This 
can be due to imperfect matching of the rf frequency 
to the spin transition. Finite Ti„h also contributes to 
matching imperfections. We have developed a theoreti- 
cal model for calculating the evolution of pab during the 



FIG. 4. Spin echo experiment. \pab\ is obtained from the 
quantity ln(J///re/)/ ln(7i/7re/), where 7i and 7/ are the trans- 
mitted probe intensities before and after the rf sequence, re- 
spectively, and Iref is a reference intensity for a known popu- 
lation difference (typically 0). The fit returns the parameters 
crA/(27r) = (2.3 ± 0.2) kHz and Tc = (172 ± 30) ^s. 



CPMG sequence. Details can be found in the Appendix. 
This model reduces to the case of the spin echo when the 
number of 7r-pulses in the sequence is set to one. The 
expression for p^b is given in Eq. p3p. The parameters 
of the model are Tc and the standard deviation of the 
fluctuation in A, cta (a global offset is also included as a 
parameter). We have fitted the model to the experimen- 
tal data. The result is represented by the line in Fig. |4l 
The agreement with the experiment is satisfactory. The 
fitting procedure yields (TA/(27r) = (2.3 ± 0.2) kHz and 
Tc = (172 ± 30) ps. Not only are these values interest- 
ing themselves, also they permit to estimate T2 even if 
the time scale of our experiment is not long enough to 
appreciate the exponential decay. The estimation is per- 
formed through the expression Tj"^ = (t'^Tc, known for 
stationary Markovian gaussian processes, and also de- 
duced in the Appendix [Eq. We get T2 = 1.01 
ms. This value doubles the one previously obtained for 
the magnetic field orientation (6* = 45°, ^ = 49°)^ The 
rise is due to the partial ZEFOZ effect discussed in the 
previous section. It is worth noting that cta is a parame- 
ter dependant on the magnetic field orientation. On the 
other hand, Tc is independent, providing us with a more 
general piece of information about the system. 

With the information gathered so far, we know that 
TT-pulses spectrally wide enough and much shorter than 
Tc are feasible in Tm^^+iYAG. This enables the design of 
CPMG sequences aimed at extending T2. In the next 
section, we will show initial evidence that this is indeed 
possible. 



B. CPMG experiments with a few 7r-pulses 

First we prepare the system by optically pumping the 
crystal into a single energy sublevel of the nuclear spin 
over an optical interval ^400-kHz broad. Then, we ap- 
ply a CPMG sequence as sketched in Fig. O^a) . It starts 
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FIG. 5. (a) Scheme of the CPMG sequence. The signes "+" 
and " that precede the label of the 7r-pulses refer to the 
pulse phase, corresponding to or tt, respectively, n and r 
stand for the number of 7r-pulses and the time interval be- 
tween them, (b, c) Coherence recovered after a CPMG se- 
quence with (b) fixed free evolution time nr = 2 ms and 
increasing n, or (c) fixed n = A and increasing r. The calcu- 
lated data are obtained from Eq. (|20[) with (7A/(27r) = 2.3 
kHz and Tc = 172 /is. 



by a 7r/2-pulse which transforms the initial population 
difference into a coherence (rotation of the Bloch vec- 
tors from the vertical axis to the transverse plane of the 
Bloch sphere). A sequence of n 7r-pulses follows, arriv- 
ing at instants tj = [2j + 1)t/2 (j = 0, . . . , n — 1) with 
phases [1 -I- {—iy^^]'K/2. They rephase the spins at in- 
stants jT, T being the time interval between two pulses. 
Just as in the spin echo experiment, another 7r/2-pulse of 
phase l)"]7r/2is applied at the end of the sequence 

for enabling the optical detection of the final coherence. 
Here again, \pab\ is recovered from the initial (right after 
optical pumping) and final transmitted probe intensity. 

We have performed two preliminary experiments to 
test the action of the CPMG sequence onto pab- In the 
first one, n was increased while the total free evolution 
time of the sequence, nr, was kept fixed at 2 ms. The 
results are exhibited in Fig. El^b). We observe that \pab\ 
rises as n is increased (t is shortened), providing evi- 
dence of the preservation of pab by the CPMG sequence. 
Our theoretical model (hollow squares) matches the ex- 
perimental data (solid squares) very well for the same Tc 
and (TA obtained from the fit of the spin echo data (Fig. 
5]). For n large, \pab\ seems to tend to a stationary be- 
havior. In fact, the rightmost data point in Fig. [Sl^b) is 
close to Tc. A plateau- like behavior is better defined in 
FiglDJc). Here, n is fixed to 4 and r is varied through 
Tc- We see that, as t is shortened, \pab\ increases until it 
reaches a plateau {\pab\ not attaining unity for t ^ is 
probably due to the reasons mentioned in the analysis of 
Fig. |4]). This behavior allows us to distinguish between 
two regimes. In fact, the borderline is in good agreement 



with the previously obtained value Tc = 172 ^s. For 
T > Tc, A fluctuates significantly during one elementary 
sequence of the CPMG process [dashed rectangle in Fig. 
[5Ja)]. In such a case, the CPMG strategy fails at keeping 
Pab- For T < Tc, the elementary sequence is short enough 
to ensure that the decoherence is suppressed or, at least, 
strongly diminished. 



C. Extension of the coherence lifetime 

Let us now turn to experiments aimed at keeping the 
coherence over long time intervals. In Fig. [B] we show 
Pab as a function of nr (total free evolution time), for 
CPMG sequences with a large number of pulses, and for 
different values of t. Each data set labeled r — 150, 100, 
10, 3 ps has been obtained by increasing n for a fixed 
value of r. In the figure we also include the standard 
spin echo measurement (labeled n = 1) for comparison. 
We observe that the data sets display a linear behavior 
in log scale. This translates the fact that the time scale 
involved is much longer than Tc [see the Appendix, es- 
pecially Eq. ((2Tt ]. The slope, defined as T2^, decreases 
significantly for decreasing t. The maximum T2 attained 
is (230 ± 30) ms, for r = 3 /is. That is over 220 times 
larger than the value obtained through the standard spin 
echo, and 450 times larger than the one measured for a 
magnetic field orientation {9 — 49°, = 45°) that does 
not satisfy any 9qA = {a = 9,(1)) condition^* (absence 
of partial ZEFOZ effect). These results prove the T2 ex- 
tension in Tm"^+:YAG provided by the application of a 
CPMG sequence with t < Tc- 
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FIG. 6. (Color onhne). Measurement of T2 through CPMG 
sequences with different r values. For data sets labeled r — 
150, 100, 10, 3 ps, T is kept fixed and n is increased. The 
data set labeled n = 1 is a standard spin echo experiment (n 
is fixed and r is increased) . The dashed lines result from least 
squares fits. 

In Fig. [3 we plot the values of T2 obtained from a 
least squares fit to the data of Fig. [5] The figure also 
displays the theoretical prediction for T2 according to 
Eq. (pij) . We see that the overall experimental and 



7 



theoretical behaviors are similar. However, the experi- 
mental T2 is always lower than the predicted one. For 
/i = 150 /iS, experimental and theoretical values are 15.4 
and 18.7 ms, respectively. However, for /i = 3 /iS, the 
corresponding values are 0.230 and 43 s. The ratio the- 
ory/experiment is not far from unity for large values of 
T, but it attains 190 for t — 3 iis. This is partly due 
to the cumulated imperfections of the spin rephasing by 
one single elementary sequence. For t ~ 150 /iS, only 100 
elementary sequences are applied to reach 15 ms of free 
evolution, but for r = 3 fj,s, 5000 sequences are necessary. 
In addition to this source of discrepancy between experi- 
ment and theory, it is worth recalling that we have made 
use of a very simple model of decoherence in virtue of 
its straightforward mathematical treatment (see the Ap- 
pendix). This model assigns decoherence to fluctuations 
of the spin transition frequency only. The latter quantity 
is defined as a stationary, Markovian, gaussian stochastic 
process, three properties that suffice for complete charac- 
terization. Spin-spin interactions between the impurity 
ion and the spin bath of the crystal (other than those of 
the flip-flop kind) can be classed into this category with- 
out introducing much error. However, spin-lattice inter- 
actions, for instance, of rarer occurrence and responsible 
for Ti (measured to be ~ 1 min), are not at all taken 
into account in our model. Flip-flop interactions are dis- 
carded as well. When looking at long times, which be- 
comes possible thanks to the DD compensation for spin- 
spin decoherence, these effects become dominant, hence 
the discrepancy with the experimental results. 
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FIG. 7. (Color online). CPMG-extended T2 as a function of 
T. The experimental values (solid black symbols) result from 
the least squares fit of the data of Fig. [B] The theoretical 
prediction (solid red line) represents Eq. (|2ip . The dashed 
green line indicates the lower bound (2Ti)~^ (Ti ~ 1 min). 
The ratio between theoretical and experimental T2 values is 
also plotted (hollow blue symbols; these data refer to the right 
vertical axis). 



IV. CONCLUSIONS 

TT-Pulse-based DD techniques on Tm'^+iYAG posed 
the problem of high rf power requirements imposed by 



the large inhomogeneous broadening of the spin transi- 
tion. We have shown that it is possible to minimize this 
broadening by carefully choosing an orientation for the 
static magnetic field such that the sensitivity of the spin 
transition frequency to the field orientation is minimal. 
This way, a large span in magnetic field orientation seen 
by spins at different sites in the lattice translate into a 
small span of transition frequencies. The inhomogeneous 
broadening has been reduced by a factor > 5. At the 
same time, T2 has been doubled as a consequence of a 
partial ZEFOZ effect. 7r-Pulses effective over the whole 
inhomogeneous line are now feasible on Tm'^+ : YAG with 
reasonable rf powers. This has enabled the application 
of a CPMG DD sequence. We have observed a clear 
increase in as the time length of the elementary se- 
quence in the CPMG protocol was reduced. A maximum 
value of (227 ± 30) ms was obtained. It represents a 445- 
fold rise with respect to previous experimental values. A 
simple theoretical model we proposed showed excellent 
agreement to our experimental data for sequences with 
few pulses (~10). It allowed the extraction of the pa- 
rameters that characterize the fluctuation: its standard 
deviation and its correlation time, estimated to (2.3±0.2) 
kHz and (172 ±30) //s (at 1.7 K), respectively. The com- 
parison between experiment and model declines as the 
number of pulses is increased. We assume this is due to 
accumulation of pulse imperfections and to decoherence 
mechanisms disregarded in the model. 
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APPENDIX 

Here we aim at calculating the evolution of the co- 
herence during a CPMG sequence. We will consider an 
inhomogeneously broadened atomic transition centered 
at Aq. Each ion in the inhomogeneous line suffers from 
crystal field fluctuations that induce fluctuations in its 
spin transition frequency. Its detuning from the line cen- 
ter varies as 

<5(t) = ^ + C(t), (6) 

where S is the time-averaged detuning. We assume the 
initial situation (t = 0) in which the coherence Pah^S- — 
0) = is the same for all the ions, and there is no pop- 
ulation difference (all Bloch vectors aligned along an axis 
contained in the transverse plane of the Bloch sphere). 
This scenario can be reached, for example, by emptying 
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one of the levels of the spin transition (optical pumping) 
and applying a rf 7r/2-pulse, as done in the experiments 
described above. 

The CPMG sequence consists of a succession of n tt- 
pulses. The 7r-pulses are applied at instants t = {2j + 



1)t/2, j = 0, . . . ,n — 1. We are interested in the echo 
intensity at t = nr. This is proportional to the coherence 
averaged over the inhomogeneous distribution, pat- As 
each TT-pulse of the sequence conjugates the coherences, 
in the reference frame rotating at angular frequency Aq 
Pab is given by 



Pab{n,T) = p 



exp 



/ /•(2"-l)f 

I I dtS{t) - [i dtS{t) - 



-i / dtSit) 
Jo 



5i / rH 

i I dtS{t) - (i dtS{t) 



(7) 



Pab{n,T) = p^^H exp<^ i 



+(-ir 



/•(2"-l)i 

dtC{t) + {-i) / dtc{t) 

(2n-l)f J(2n-3)f 



dtQit) + (-1)" / dtQit) 
Jo 



+ --- + {-iy 



dtc{t) 



(8) 



where g{S) stands for the inhomogeneous distribution 
and (...) represents the statistical average. For going 
from Eq. ([7]) to Eq. ([5]) we have taken into account that 
the contribution from 5 cancels out. Equation ([Sj can be 
written as 



with 



(9) 



r 



With this in mind, the expression for the coherence re- 
duces to 



Pabin,T) = p'-Je 



- JO),-7(",^) 



with 



7(n,r) 



dt'sr{t)C{t) 



(13) 



(14) 



(—1)" < i < -, 
Mt) = { (-!)"-■'■ (2j-l){'<t<(2j + l)§ , (10) 
1 (2n-l)§<^<nr 

where j = l,...,n — 1. 

Let C,{t) be a stationary, gaussian, Markovian pro- 
cess. According to Doob's theorem)^ only the Ornstein- 
Uhlenbeck'^* process simultaneously satisfies those three 
assumptions. The corresponding autocorrelation func- 
tion reads 



(C(t)C(i')) = ^Ae-l*-*'l/^% 



(11) 



Let us perform the calculation for 7(n, r). 



riT /•nr 



j{n,r)^-l dt I dt'sr{t)sr{t'mt)at')), 



fo Jo 

2 pnr pnr 



nnr nnr 

dt dt'sr{t)sr{t')e 
2 Jo Jo 



'^--l*-*'l/^%(15) 
(16) 



The above double integral simplifies if we consider the 
relation Jldt'Jldt" fit', n = 2 Jljt' dt" f{t' X), 
valid for any f(t',t") symmetric after permutation of t' 
and t" . We therefore get 



where a a and Tc represent the standard deviation of the 
fluctuation and the correlation time, respectively. Since 
C(i) is a gaussian process, it can be shown that 



j{n,T)^al dt dt'sr{t)sr{t')e-^'-*">/^', {17) 
Jo Jo 



= a\ I dt Srit)hr{t), 

Jo 



(18) 



with 
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Kit) = i-inaATcVe-- { 



e- - 1) 
- 1 



g2xc - 1 



e + 1 L 

^"-1 ^ ^l + (-l)"e("-i)^ 
(2i-i)2^ 



e^c + 1 

+ (-!)"(' 



e^c — e' 



0<t< 



(2j-l)§ <t < (2j + l) 



(2n- l)f <t<nT 



(19) 



where j runs from 1 through n — 1. By computing the 
integral in Eq. (ITSt . we obtain the final expression for 
7(n, r): 



1 2 / T 
tanh — 




1 — sech 

2Tr 



(20) 



It is understood in Eq. ([^0]) that t = nr. 

We see from Eq. ([^0]) that the coherence does not de- 
cay in an exponential manner. However, the exponential 
behavior is recovered in the i 3> Tc limit with a charac- 
teristic time given by 



T2 ^(t) = cr^Tc 



r 



2rc 

tanh , 

T V 2t, 



(21) 



Moreover, if r ^ Tc, 12 reduces to 



1 ^2 ^2 

„ 1 / s t (T A T 

LZ Tr 



(22) 



From Eq. ((20|) . we can easily calculate the coherence 
lifetime obtained by a single spin echo by assigning n = 1 
and T ^ t. This yields 



7.e(i) = (fTATe)' (^^ + Ae'^ - e"^ - 3 ) . (23) 
In the limit t ^ Tc we get the standard result 



(24) 



It is worth mentioning that the development described 
here to obtain 7(71, r) is analog to the extension to n tt- 
pulses of the calculation performed by Herzog and Hahn 
in Ref. [s^ for one 7r-pulse. The same result is found. 
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